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Motivated by recent discovery of ferromagnetism in cubic perovskite BaFeOa under small magnetic 
field, we investigate spin order in BaFeOa and isostructual SrFeOa by the first principles calculation. 
The on-site Coulomb and exchange interactions are necessary for the helical spin order consistent 
with experiments. SrFeOa exhibits stable G-type helical order, while A- and G-type helical orders in 
BaFeOa are almost degenerate at short propagating vector with tiny energetic barrier with respect 
to ferromagnetic spin order, explaining ferromagnetism under small field. The results are consistent 
with model calculation where negative charge-transfer energy is explicitly taken into account. 

PACS numbers: 75.30.-m, 75.30. Et, 75.50.Bb, 75.40. Mg 
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I. INTRODUCTION 

The coexistence of spin order and metallic conductivity 
in Sd transition-metal oxides have attracted a lot of re- 
search attention since the discovery of high-temperature 
superconductivity in copper oxides. The metallic con- 
duction is accessible by the introduction of carriers into 
either Mott-Hubbard or charge-transfer insulators in the 
Zaanen-Sawatzky- Allen phase diagramii Without carrier 
doping, metallic conduction is expected in the phase dia- 
gram when systems are located in the area with small or 
negative charge-transfer energy A. Such a situation can 
be obtained for large atomic number and high valence of 
id transition metal. 

A typical example of the negative A compounds is cu- 
bic perovskite SrFcOs, where formal valence of iron is 
Fe'*"'" (3d^) and the effective value of A is estimated to 
be A '^ — 3 eV<^ Metallic conductivity has been pre- 
served even below an antiferromagnetic (AFM) transi- 
tion temperature Tn -- 134 K,^"^ Below Tn, SrFeOa 
shows a G-type helical spin order whose propagation vec- 
tor is parallel to the [111] direction. The observed vec- 
tor is q = 0.112(1,1,1) X 27r/as with lattice parameter 
as = 3.85Ai^ The magnetic moment is 3.1/iB per iron 
at 4.2 Kii Recently versatile helimagnetic phase diagram 
under magnetic field has been established and unconven- 
tional Hall resistivity has been discussed in connection 
with topological Hall effect in noncoplanar spin texture 
with scalar spin chirality;^ 

In contrast to SrFeOa, CaFeOs is insulating below 
290 K)^ where a charge disproportionation by the for- 
mation of Fe^"*" and Fe^+ emergesJ^ The difference may 
come from the compression of lattice due to small ironic 
radius of Ca^"*" compared to Sr^+. It is thus interesting 
to examine oppositely the effect of lattice expansion by 
replacing Sr^+ by Ba^"*". 

Very recently cubic BaFcOs has been synthesized by 
low-temperature chemistryjii The ground state seems to 



be AFM and there is no charge disproportionation and 
no structural distortion down to 8 K, similar to SrFeOs. 
However, its helical order is different from SrFcOa: A- 
type helical order with a propagation vector along the 
[100] direction. Here, q = (0.06, 0, 0) x 27r/aB with lattice 
parameter ae = 3.97A. By applying small magnetic field 
~ 0.3 T, the AFM state changes to ferromagnetic (FM) 
one whose saturated moment is 3.5 /^e/Fe. This is also 
different from SrFeOa where saturated ferromagnetism 
(3.5 /^b/Fc) is achieved by applying more than 40 Tj^ 

The helical spin order in iron perovskites has been ex- 
amined theoretically based on a double exchange model 
explicitly including oxygen 2p orbitalj^ By adding su- 
perexchange interaction Jse between localized spins, the 
G-type helical order is stabilized. Though the reduction 
of Jse changes the G type to the A type in small and neg- 
ative A region, it is not clear whether the change really 
corresponds to the difference of real materials between 
SrFeOs and BaFeOs. In order to include material infor- 
mation in theory, we need to perform the first principles 
calculation on magnetic structures in these compounds 
based on density functional theory (DFT). 

In this paper, we perform non-coUincar spin polarized 
DFT calculations in SrFeOa and BaFeOs as a function 
of the propagation vector q for helical spin order. We 
find that the DFT calculation within local spin density 
approximation (LSDA) requires on-site Coulomb interac- 
tion U and exchange interaction J to explain experimen- 
tally observed propagation vectors. The G-type helical 
order in SrFcOa is nicely reproduced by the LSDA-f t/ 
with L/ = 3 cV, J == 0.6 cV. In BaFeOs, although the 
G-type order is lower in energy, the energy difference 
between A- and G-type orders gets extraordinary small 
around experimentally observed propagating vector in 
contrast to SrFeOs. The energy difference between fi- 
nite q and q = (FM) orders is larger in SrFcOs than in 
BaFcOs, consistent with experimental observation that 
SrFcOs requires a larger magnetic field to achieve satu- 



rated FM state. The difference of the two compounds 
is attributed to tfie difference of lattice constant. This 
is confirmed by performing a model calculation that in- 
cludes both negative A and the difference of the lattice 
constant through the distance dependence of hopping and 
superexchang parameters. We also find that the density 
of states (DOS) just below the Fermi level is dominated 
by oxygen 2p component, consistent with the view of 
negative A. 

This paper is organized as follows. In Sec. II, we in- 
troduce the computational method for helical spin order 
using DFT. The results by DFT calculation are shown 
in Sec. III. In Sec. IV, we introduce a double exchange 
model, and show results consistent with the DFT calcu- 
lations. We summarize our results from the DFT and 
model calculations in Sec. V. 



II. COMPUTATIONAL METHOD 

In order to treat helical spin order, we must constrain 
the form of wave function J^ii^ The wave function ^k(r) 
under the propagation vector q could be written as 
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*k(r) 



gj(k-q/2)T 







i(k+q/2)-r 



/4(r) 



(1) 



where /i^(r) is translational invariant Bloch wave func- 
tion with spin a before transformation. To determine 
the optimal propagation vector q, we calculate the q de- 
pendent total energy E{q) per unit cell by VASP code J^ 
Non-collinear spin polarized calculations are performed 
within LSDA with and without U and J.— The calcu- 
lation is done in primitive cell with 10x10x10 fc-points 
and energy resolution 0.01 meV per cell. Since the or- 
bital moment is quenched completely by our calculation, 
the spin-orbital coupling effect is not under considera- 
tion. The cubic crystal structure (space group Pm-3m) is 
used. Crystal distortions may happen resulting from the 
spin lattice coupling, but this distortion is not detected 
by experiment down to 8 K. Wc thus consider that the 
distortion is very small and ignore it in our calculations. 



III. RESULTS BY DFT CALCULATION 

Figure [T] shows the total energy of SrFeOs as a func- 
tion of the propagation vector by LSDA, where (j) in the 
horizontal axis is defined by q = 0(1, 1, 1) x 27r/as for 
the G type and q = 0(1,0,0) x 27r/as for the A type 
with as = 3.85Ai^ The total energy is measured from 
the = (FM) state: AE{(j)) = E{<j)) - E{(j} = 0). 
Within LSDA, the ground state is the FM state. If 
we take C/ = 3 eV and J = 0.6 eV, the minimum en- 
ergy of SrFeOa shown in Fig. 2 is located at = 0.11, 
which agrees with the experimental findings. The mag- 
netic field required to overcome the energetic barrier of 




FIG. 1. (Color online) The dependence of total energy 
difference per unit cell, AE{4>) = E{4>) — E((j) — 0), obtained 
by LSDA in SrFeOa. For the definition of 0, see text. 



A£'(0 = 0.11)=:3.0 meV is approximately equal to 52 T, 
which is close to the experimental magnetic field (^ 42 T) 
at which the FM moment is saturated)^ The calculated 
moment at = 0.11 inside the muffin-tin spheres is 
3.0/iB/Fe, close to the experimental value of 3.1//B/Fe.— 
A£:(0) for BaFeOs with lattice constantii ob = 3.97A 
is shown in Fig.[2l where the same U and J are employed. 
We find that the total energy is almost independent of 
wave vector in both A- and G-type helical spin orders 
when < 0.08. This is consistent with the fact that the 
magnetic field to obtain saturated FM state is smaller— 
(~1 T) than that in SrFeOa. The energy difference be- 
tween A- and G-type is tiny, so the A-type helical spin 
order may be stabilized by, for example, introducing cor- 
relation effect that is beyond the present calculation. The 
small energy difference may affect physical properties at 
finite temperatures. 

In order to clarify the electronic structure of BaFeOa, 
wc show in Fig. H the DOS for the FM state (0 = 0), 
since this state is very close to the stable helical state. 
Negative DOS represents down-spin component. The 2p 
electrons of oxygen and 3d of iron hybridize intensely in 
energy scale —4 '^ eV, and spin-up component crosses 
the Fermi level, which indicates BaFeOs is half-metallic. 
The fact that DOS near the Fermi level is dominated by 
oxygen 2p is consistent with the view of small or negative 
A in high valance 3d transition-metal oxides. We also 
find that the i2g orbitals arc almost occupied by up-spin 
electrons but empty for down spins. This indicates a 
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FIG. 2. (Color online) The (j) dependence of total energy 
difference per unit cell, AE{cj>) = E{4>) — E{(j) — 0), obtained 
by LSDA+t/ in SrFeOs and BaFeOa. For the definition of 4>, 
see text. 



localized nature of t2g spins. 

What is the origin of different stable helical structures 
in SrFeOa and BaFeOa? The most significant effect may 
come from lattice parameter from the results of our DFT 
calcination. Since ionic radius of Ba^"*" is larger than 
Sr^+, BaFcOa can be considered as expanded SrFeOa. 
We have confirmed that the energy difference between 
G type and A type in SrFeOa decreases with increasing 
the lattice parameter (not shown). To get better under- 
standing on the helical spin order and the effect of lattice 
in SrFeOa and BaFeOa, we connect our DFT calculation 
with a double exchange model in Sec. IV. 



IV. RESULTS BY MODEL CALCULATION 

Mostvojii^ has shown that a double exchange model 
including explicitly oxygen 2p orbitals bears the A-type 
helical spin order of localized t2g spins under small or 
negative A. By adding superexchange interaction Jse 
between the localized t2g spins, the G-typc helical order 
emerges. The Hubbard U in our first principle calcula- 
tions is crucial for stabilizing the t2g localized spins so as 
to be consistent with the double exchange model. This is 
realized in Fig. [3]where all of up-spin electrons of the i2g 
are deep in energy. In addition, we expect superexchange 
process by making use of U, though the precise estimate 




Energy(ey) 

FIG. 3. (Color onhne) Density of states (DOS) of BaFeOs in 
the FM state calculated by LSDA-f f/ with U = 3.0 eV and 
J = 0.6 eV. The positive side of DOS denotes the up-spin 
DOS, while the negative side denotes the down-spin DOS. 
Dotted vertical line at zero energy represents the Fermi energy 
Ef. (a) DOS of Fe 3d t2g and eg orbitals. (b) DOS of O 2p 
orbitals. 



of the process is difficult. We note that small or negative 
A is necessary for the double exchange model to produce 
both the G- and A-type orders. In our LSDA+U calcu- 
lations, the presence of 02p DOS just below the Fermi 
level (see Fig. |31[b)) will be an indirect evidence of small 
or negative A. 

In order to make clear the correspondence between the 
double exchange model with Jse and our LSAD+U , we 
have to show, at least, that the double exchange model 
can explain the difference in SrFeOs and BaFeOs. For 
this purpose, we examine the model following Mostvoy's 
procedureJ^ By assuming that the localized t2g spin at 
site i is in the y-z plane, i.e.. Si = 5'(ysinq-Xi -|-zcosq- 
Xi) with unit vectors y and z along the y and z directions, 
respectively, the effective Hamiltonian of the model may 
be written aai^ 



k,a,(5 k,(5.(T 

+JsE /_^ Si • Sj, 



(2) 



with 



PkSa = 2 



qs ks . qs . ka 

cos — COS —Pksa - sm — sm —pks- 



where a runs two €„ orbitals of 3z 



(3) 
■^ and x^—y'^, S runs 
three component of x, y, and z, (i,j) runs the nearest- 
neighbor pairs, dka represents the annihilation of Cg elec- 
tron, and PkSa represents the annihilation of 2ps electron 
with spin a. The hopping matrix clement ta^s is given 

tx^-y^,x = —tx2-y2^y = v3(pdo')/2, and tx2_y2., = 0, 
with the hopping parameter (pda). We ignore oxygen- 
oxygen hopping for simplicity. We find from ^ that 
finite q mixes both down and up spins for holes on oxy- 
gen. 

We calculate the total energy of the Hamiltonian ([2]) 
with one eg holei^ and plot it as a function of 4> for both 
the G- and A-type helical orders as shown in Fig. S) For 
SrFeOs, we take A = —3 eV and {pda) = 1.3 eV ac- 
cording to photoemission analysis.— We determine a re- 
maining unknown parameter JgE in order for the energy 
minimum to be located at the experimentally observed 
propagation vector, as shown in Fig. |4l^a). Jse is, then, 
estimated to be Jse = 9.1 meV. Using the same param- 
eter set, we find that the minimum energy of the A- type 
order is higher by 0.92 meV and the energy of the FM 
state {(j) = 0) is also higher by 1.66 meV. The latter en- 
ergy corresponds to the magnetic field of 28.7 T that is 
slightly smaller than the FM transition field ~42 T for 
SrFeOg-^ 

For BaFcOa, we use a slightly smaller (pda) = 1.17 eV 
taking the lattice expansion into account. By minimizing 
the energy at around experimentally observed (j), we find 
Jse = 7.34 meV, which is smaller than that for SrFeOs 
and is a reasonable value as a result of the expansion of 
the lattice constant compared with SrFeOa. The energy 
relative to (/> = is very small with 0.017 meV as shown 
in Fig. mjb), which corresponds to 0.3 T that is compa- 
rable to the experimental value of 1 T. Such a small en- 
ergy difference mainly comes from the decrease of {pda) 
as compared with SrFeOs. We note that the reduction 
of JgE makes the energy difference between the G- and 
A-type orders considerably small, and the difference be- 
comes almost zero for BaFeOa. 

We realize that the results of the model calculation 
are qualitatively consistent with the LSDA-|-[/ results. 
Quantitative differences may be due to many assump- 
tions for the construction of the double exchange model. 
Nevertheless, the agreement between the model calcula- 
tion and LSDA-|-[/ is satisfactory. This means that the 
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FIG. 4. (Color online) The dependence of the total energy 
difference relative to the = 0, AE{(j)), of the model ((2}. 
(a) SrFeOa and (b) BaFeOa. Parameters denoted in each 
panel are taken in order for the minimum of the energy to be 
located at experimentally observed <j) (see text). The G(A)- 
type helical spin order is represented by solid black (broken 
red) line. 



LSDA-|-f/ results contain the essential part of the model 
calculations, i.e., negative A and competition between 
the double exchange and superexchange processes. 



V. SUMMARY 

Wc have examined helical spin order in cubic pcr- 
ovskite SrFcOs and BaFeOs by the first principles cal- 
culations based on DFT. We have found that U and J is 
necessary for explaining experimentally observed prop- 
agation vectors of spin. Including U and J to LSDA, 
we have obtained the G-typc helical order for SrFeOa 
and almost degenerate A-type and G-type spin order for 
BaFeOa- The energy difference between the minimum- 
energy state and FM state is larger in SrFeOa than in 
BaFeOs, consistent with experimental observation that 
SrFeOs requires larger magnetic field to achieve satu- 
rated FM state. The difference of the two compounds is 
attributed to the difference of lattice constant. We have 
shown clear correspondence between the first principles 
calculation based on LSAD+L/ and the double exchange 
inodel implemented by superexchange interaction, con- 
firming the importance of the competitions between dou- 
ble exchange and superexchange interactions. The cor- 
respondence also implies that the characteristics of neg- 
ative A is included in the first principles calculations. 
Correspondingly, the DOS of oxygen 2p component in 
FM BaFeOs is just below the Fermi level, which is con- 



sistcnt with the view of negative A. 
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